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Abstract 

We present a new method for constructing solutions to nonlinear evolu- 
tionary equations describing the propagation and interaction of nonlinear 
waves. 

In the present paper, with the help of some simple examples, I demonstrate 
a new approach to the construction of asymptotic solutions to differential equa- 
tions. I choose very simple examples (the Hopf equation and its multidimen- 
sional analog) and, on purpose, omit details in the proofs of estimates (which 
are obvious here). 

Usually, by saying that a function is an asymptotic (approximate) solution 
of a differential equation, we mean that this function satisfies the equation with 
a small discrepancy. The smallness of the discrepancy is understood as the 
smallness in some uniform metric under the assumption that a small parameter 
tends to zero. 

A function is called a weak asymptotic solution if, after the substitution 
of this function into the equation, there is a discrepancy that is small in the 
weak sense as a small parameter tends to zero. In this case the functionals are 
assumed to depend on time as on a parameter. 

For example, under this approach, the ^^-approximation of a generalized 
function turns out to be its weak asymptotics and we can choose generalized 
functions to be the initial conditions and use their approximations for construct- 
ing the solutions. In this case, we obtain a small parameter, which is either the 
parameter of approximation or the small parameter in the original equation. In 
the latter case, this original small parameter is taken to be the parameter of 
approximation. 

In fact, this approach is close to the ideas proposed by J. F. Colombeau and 
other authors who constructed different algebras of generalized functions. The 
method itself was first introduced in our papers with V. M. Shelkovich initiated 
by the works of J. F. Colombeau and M. Oberguggenberger with coautors and by 
discussions with J. A. Marti, V. V. Zharinov, and S. Pilipovich. The difference 
with the traditional approach is that in our approach the mollifier is chosen not 
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from the consideration of the algebraic construction but from the consideration 
of the original differential equation. 

In some cases (shock waves), the solution is independent of the choice of 
the mollifier, while in other cases (solitons, kinks) the solution depends on this 
choice. 

If the original equation contains a small parameter, then we, in fact, deal 
with rcgularizations by small viscosity or small dispersion. In this case, to 
calculate a weak asymptotics, we need to calculate the zero viscosity or zero 
dispersion limits. Hence we arrive at the problem of constructing a definition 
of a weak solution which admits this passage to the limit. 

It turned out that the approach developed here can be used for describing 
both the propagation of nonlinear waves and, which is the most important, 
their interaction. It is well known that the problem of interaction of nonlinear 
traveling waves (for instance, of two kinks or two solitons) in the case of a 
single spatial variable can be formulated as a problem of constructing the exact 
solution of a nonlinear equation with two spatial variables corresponding to the 
coordinates of the wave fronts. If the initial problem can be integrated by the 
method of the inverse scattering problem or in any other way, then one can 
write out the solution of the above-mentioned equation explicitly, which allows 
one to describe the interaction analytically. In our approach, to describe the 
interaction, one needs to solve an ordinary differential equation (or a system 
of such equations) with a small parameter. Solutions of such equations can 
be constructed by using a well-known technique. In what follows, we consider 
the main technical tools and some examples which allow us to demonstrate the 
abilities of our approach. 



1 Some weak asymptotic formulas 

(a) Let ijj(z) € §(R 1 ), where S is the Schwartz space. We consider the function 
-w{{x — a)/s) and calculate its weak asymptotics. Treating -u)((x — a)/e) as a 
generalized function, for any function rj{x) G C§° we have 

— -^j,r](x)) = — J — -^jr](x)dx = J w(z)r](a + ez) dz 

fc>0 

where the last relation is formal and means that the left-hand side can be 
represented as the asymptotic series given on the right-hand side, 

O fe = J iv(z)z k dz. 

We define by Ox»(s a ) and element of V such that 

f(x,s) = O v ,{e a ) & (f(x,e), V (x)) = 0(e a ), 
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where the last O-estimate (which must hold for any function r)(x) G Cq°) is 
understood in the usual sense. Then for any N we can write 



1 fx — a\ x „ e 

• : - k \ 

k=0 



N k 

n k £ -{-l) k 6^{x -a) + O v ,(e N+1 ). 



(b) Let u)i(z) , u>2(z) G ^(K 1 ). Let us consider the weak asymptotics of the 
product u)\((x — ai)/e)u)2((x — a 2 )/s). We have 

fx-a 1 \ (x-a,2\ , f (x-a\\ (x-d2\ , . , 

— - — — - — ),V{x)) = / — - — — - — Ji]{x)dx 

= srf(ai) J uji{z)ijJ2(z dz + 0(e 2 ) 

= erj(a 2 ) J Wi {z + lo 2 (z) dz + 0(e 2 ), 
Aa = CL2 — a\- 

Finally, we obtain the following formula that is uniform and symmetric in a\, 0,2'- 
^{^iy 2 (^^) = \[e5{x 0l ) + sS(x 02)]S(^) + 0^(e 2 ), (1) 



where 



b(^——^J= J ui{z)lu2{z — ^ dz = J uji(^z + —^jLo 2 (z) dz. 

(c) Now let wi(«),W2(«) G C°°, ^ G S0R 1 ), lim^.^Wi = 0, lim^^u;; = 
1, i = l,2. 

Let us calculate the weak asymptotics of the derivative 

1. fx — a\\ fx — ci2\ 1 fx — ai \ . fx — a 2 \ 
= 7 Wl (-i-j wa (-i-j + i Wl (-i-j wa (-i-j- 

Just as previously, we have 

1. /x — ai\ fx — a 2 \ 1 /a; — ai\. fx — a 2 \ 

= 5(a:-ai)Bi(^) + ^ - a 2 )B 2 (^) +0 D /(e), 



where 



Bi(^^ = j ui(z)u2(z - dz, B 2 (PP) = J ui(z+ ^P)w2(z)dz, 
Bi(oo)=0, Bi(-oo) = l, Bi(z) + B 2 («) = 1. 



We have 
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Calculating the primitive, we obtain 

Aa\ „, N „ / Aer 



. 1 (^y 2 (^)=6( X -a 1 )B 1 (^) + e( X -a 2 )B 2 (^) + OMe). (2) 



(d) Under the assumptions of item (b) and the condition that J Wi(z) dz = 
1, the functions Wi((x — a^/e) are approximations (weak asymptotics) of the 
functions ed(x — di), 

Ui((x - a l )/e) = e5 £ ^(x - a;) 
Hence we can rewrite (1) as 

eS £i i(x ~ ai)e5 £ ^{x - a 2 ) = ^[eS(x - a x ) + sS(x - a 2 )]B(Aa/e) + O v >(e 2 ). 

In a similar way, under the assumptions of item (c), u>i((x — a\) / e) = 9 £ .i{x — ai) 
are approximations of the Heaviside ^-function. Hence we can rewrite (2) as 

9 e .x{x - a 1 )9 ea {x - a 2 ) = 0(x - a 1 )B 1 {Aa/e) + 9(x - a 2 )B 2 (Aa/e) + C /(e). 

2 Nonlinear structures 

We show how the above formulas can be used to describe interaction of nonlinear 
structures. 

(a) Interaction of shock waves for the Hopf equation. Let us consider the 
Cauchy problem 

L[u] = u t + {u 2 ) x = 0, u| t=Q = M + Ui9(-x + oi) + u 2 6(-x + a 2 ), 

where Ui are positive constants, a 2 < a\. We approximate the initial condition 
according to the formulas from item 1 (c) and seek the weak asymptotics of the 
solution in the form 

u e (x,t) = u + u 1 e ,i{-x + ipi(t,s)) + u 2 9 £ ^ 2 (-x + ip 2 (t,e)), 

(£i(0) = ai, <^2(0)=a 2 . 

Calculating the weak asymptotics of the expression (u e ) 2 according to the for- 
mulas from item 1 (c), we obtain 

(u e ) 2 = ul + (ul + 2u u 1 )9(-x + cpt) + ul + (uj + 2u u 2 )9(-x + ip 2 ) (3) 



+2u\u 2 
where 



-x + i Pl )B 1 + 9(-x + V2 )B 2 (-^) 



Bi(A(p/e) = f oji(z)lj 2 (z + Aip/e) dz, 

B 2 {Aip/e) = J cji(z - Aifi/e)d) 2 (z)dz, Atp = ip 2 - ipi, 
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and, in contrast to item 1 (c), we have oo) = 0, £>i(oo) = 1, but as before, 
B!+B 2 = l. 

We substitute the approximation of u E (x,t) into the equation and require 
that the relation L[u s ] — Oxv(e) must be satisfied (this is the definition of the 
weak asymptotics solution in this case). Moreover, the function L[u e ] must be 
weakly piecewise continuous with respect to t for each fixed e. We obtain 



d A 
L[u £ ] = ^2 [ Uk ~d~f ~ 2u ° ?ife ~ u t- 2u 1 u 2 B k ^-^j S(-x + ipk) + V '{e) 



-r,-—— ~ 2ll,n1I,u — ?l? — 2ll,i1l.tRu( — —\ 

k=l 

Hence, in view of the definition of the weak solution, we have 



~dt~ 



2u + u k + 2u 3 - k B k (—^y k = l,2. (4) 



For Aip < (before the interaction) we have Bi(A<p/e) — and B 2 {Aip/e) = 
1 up to 0(e N ) and Eqs. (4) describe the propagation of noninteracting shock 
waves. We write 

ifW = (2u + U\)t + Ol, (f20 = {2(UQ + Ui) + U 2 )t + 02, 

then ipo(t) — (P20 ~ Pio is the distance between the fronts of noninteracting 
waves. At time t* , ipo(t*) = 0, the fronts merge. To construct a formula that is 
uniform in t and represents a weak asymptotic solution, we shall seek the phases 
(p k (t, e) of shock waves in the form 

<Ph{t,e) = ifko(t) +V>o</>fci (t), (5) 
where r = ipo (t) je and it is assumed that 



0, fl 
Or 



= (r- 1 } 

|r | — >oc 



Calculating the limit values of <frki(°°) = (p^, obtain formulas that describe the 
coordinates of the fronts of shock waves <p\ (t) after the interaction. Substituting 
expressions (5) into Eq. (4), we obtain 

~dT ~A^ TLf> ^ = 2u + Ul + 2u 3 - k B k ^— J, fc= 1,2. 
We calculate the difference of these equation: 

¥=F(p), p=^, F(p) = 2B 2 (p)-l. 
dr s 

The boundary condition for this equation has the form p/ T \ r ^_ ryo —* 1- The 
equation F(p) = has a single root po and B2(po) = £?i(po) = 1/2, which 
implies that after the interaction (ipo > 0, r = tpo/s — > 00) the wave fronts 
move with the same velocity 

dwt 

-^f = 2u + 1*1 +it 2 , A: =1,2. 
at 
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Figure 1: 

From (4) for the functions 4>ki we obtain 

6n = (-i)^ x 7 2 T\ [\B 2 {p)-i]dT>. 

The weak limit Uo(x, t) of the weak asymptotic solution u e (x,t) satisfies the 
classical definition of the generalized solution (in the form of integral identity) 
and the stability condition. 

(b) Interaction of weak discontinuities. Generation and decay of shock 
waves. 

We again consider the Hopf equation and pose the following initial condition 

U lt=0 = U + u l( a l ~ x )+ - "l( a 2 - x) + , 

where a\ > a 2 , z + = z8(z), = const > (see Fig. 1). 
We shall seek the weak asymptotic solution in the form 

u E (x,t) — u + ui{ipx{t, e) - x)9 £ s{— x + (fi(t, e)) 
-u 2 (ip 2 (t, s) - x)6 £t2 (-x + ip 2 (t, e)). 

In this case, the equations for the functions Ui = Ui(t,s) and ipi = ipi(t,s) are 
derived by using a somewhat different technique than that used for studying 
shock waves. Substituting the approximation of u £ (x,t) into the equation and 
taking into account the definition, we obtain 

(ux(<pi - x)+) t - (u 2 {<p 2 - x)+) t + (i4((p! - x)\) x + {u\{^ 2 - x)\) x 
+2[u wi(<pi - x)+] x - 2[u u 2 (ip 2 - x)+] x 
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-2[uiu 2 (vJi - x)(tp 2 - x)6(<pi - x)] x B 1 (Aip/e) 

-2[uitt2(vi - z)(y> 2 - a:)^(^2 - a;)] x B 2 (Aip/e) = 0, = p 2 - ipi- 

Let us consider the domain ip 2 < x < f\. We obtain 

u lt (<pi -x) +unpu + 2[u u 1 (ip 1 -x)] x - [ujitpx ~xf] x 
+2uiu 2 {(fi - x)B 1 + 2uiu 2 (ip2 - x)B 2 = 

We set x = ifi. Since we have = in our example, we obtain 

+U\ip\t — 2uqUi + 2u\u 2 AipB\ = 0. 

Substituting this relation into the last equation, we arrive at the following equa- 
tion for the function u\: 

uu — 2u\ + Auiu 2 Bi = 0. 

In a similar way, considering the domain — oo < x < tp 2 , we obtain the other 
two equations 

ip 2t -2u + 2u 1 AipB 2 (Aip/e) = 0, 
u 2t + 2u\ - 4uiu 2 B 2 (A(f/e) = 0, Aip = ip 2 — ipi. 

Let Alp < 0, then, up to 0(e N ), we have Bi(A(p/ e) — 0, B 2 (Aip/e) — 1 and 
obtain the following system of equations describing the evolution of the broken 
line until it turns over: 

(ipioYt - 2u Q = 0, (P2o)'t - 2«o + 2u w (ip 2Q - ip 10 ) = 0, (6) 
(<Pw)' t - 2(mi ) 2 = 0, (u 20 )' t + 2ul - 4uioM 2 o = 0, 

Solutions of this system have the form 

uio(t) = «20(*) = «?/(! - 2tu?), 
<f 10 =Oi + 2u Q t, (f 2Q = a 2 + 2[u?(ai - a 2 ) + u ]t. 

We write ipo = ip 2 o{t) — fio(t)- At time t = t* such that Vo(**) = the weak 
discontinuities merge and a shock wave is generated. 

To construct formulas that are uniform in t and describe the confluence of 
weak discontinuities and the generation of a shock wave, we seek the solution 
of Eqs. (6) in the form 

<Pk{t,e) = <pko(t)+ipo</>ki{r), r = Vo/£, 

where the functions <f>ki (r) = satisfy the same conditions as in item 2 (a) 
We shall seek the functions u £ (t, e) in the form 

u k {t,e) = il>o{0)ui/(^o + eg k (T)). 
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Here we assume that the functions <?fc(V) behave in the same way as the functions 
4>ki{T) and take into account the relation 

uio/u? = 1/(1 - 2tu?) - MWMt), 

follows from the equation tpo t + 2^10^0 = 0. After simple calculations we see 
that the function g = g\ = g 2 satisfies the equation g + 2(1 — B 2 (p)) = and 
the function p = p(r) = Atp/e is a solution of the boundary problem 

p=l-2B 1 (p), p/r^^-l. 

As before, the equation 1 — 2B\(p) has a single root p = p such that Bi(p ) = 
B 2 (po) = 1/2 and p — > po as r — > 00. This allows us to calculate the solution 
for Aipo > (i.e., after the interaction) or as r — > 00. 

We introduce the function G(t) = T+g(r). Obviously, G — p, G/t| t _ > _ oo — ► 
+ 1, and we choose 



/oo 
(l-2Bi(p))dT / +p . 
-00 



On the other hand, we can express the functions Ui via the function G: 

eG spo 

We calculate the limit {fk)t as t ^ 00 of the velocities of the weak discontinu- 
ities 

, a+ n 2</>o(0)u? 1 „ . n 

((^ 2 ) t = 2 «o ~£Po = 2w + (ai - a 2 )u^, 

, x+ „ 2^o(0)u? 1 N n 

(Vi)t = 2"o 7£po = 2u + (ai - o 2 )ui, 

epo 2 

which coincides with the velocity of the shock wave 

U(x,i) = u + (en - a 2 )u?0(-a; + </? + (i)), 

where (p + = (f 2 )t = (ViT)t- By using the explicit formula for the solution 
u e (x,t), we can easily show that 

w — lim u e (x, t) — U(x, t), t>t*. 

e— >0 

To this end, we rewrite the above-constructed solution u s (x,t) in the form 

U e (x,t) = U + U 1 (f 1 -if 2 )de,l(fl + Ui(x-<P2) [^£,2(^2 ~X) -0 e ,l(Vl -X)]. 

Consider the second term. We have 

«i(¥>i - ^2) = ^ = -0o(O)ui = (01 - a 2 )wi = U . 
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Figure 2: 



Since <^i| t>t , ~ <y9 + , the first two terms pass into the shock wave U(x,t) for 
t > t* . Consider the last term 

Ul(x-<P2)[0 e ,2(<P2-x)-0 e ,l(<Pl-x)] = Ul(x-lf 2 ) 

As was already shown, the coefficient of the expression in braces is a constant. 
The expression in square brackets is an approximation of the 5-function at the 
point (f2- Hence the entire expression in braces is small (in a uniform metric) 
as e — > 0. 

We study the problem in which a shock wave is generated by a special (piece- 
wise linear) initial condition. The case of a general smooth initial functions can 
be treated similarly. Here we need to consider a family of linear interpolations of 
this initial condition and to use the above technique on segments of the broken 
line. 

To study this problem in more detail, we note that we have considered only 
one possibility of evolution of the broken line, namely, formation of a step. 
Another mechanism of evolution is as follows: segments of the broken line are 
added to the step that has already been formed. This is the confluence of a 
weak discontinuity and a shock wave. 

Now we again consider the Hopf equation in order to study this mechanism. 
The initial condition corresponding to this type of interaction has the form 

u\ t _ = w[5#(a5 — x) + u1(ai — x)9(ai — x) — uj(a2 — x)9(a2 — x), 

where Uq, u\ are positive constants and a\ > a 2 (see Fig. 2). 



0egVP2 -x)- £ ,i(y>i - x) 



9 



Just as before, we construct the approximation of the solution in the form 

u e (x,t) = u 6 eA ((pi - x ) + «i(<Pi - x)0e,i(<Pi - x)- u x (y 2 - x)e E ^{f2 - x), 

where Ui = Ui(t,e), <pi = ipi(t,e). Substituting this expression into the Hopf 
equation, we obtain the system of equations (cf. above) 



<Pit ~ mo + 2u 1 ^pB 1 = 0, 

u lt - 2u\ + A:u\u 2 B x = 0, 

if2t - 2u B 2 + 2uiipB 2 = 0, 

u t - u ui(l - 2Si) = 0, 



(7) 



where Bi = Bi(Aip/e) are the functions derived above, Aip = ip 2 — ipi. Before 
the interaction, we have ip 2 < ipi, Aip/e ~ — oo, and B\ = 0, B 2 = 1 with 
arbitrary accuracy in e. Denoting by ip w , u w , <p 2a , u 00 the solution of system 
(7) with Bi = 0, B 2 = 1, we obtain the following system of equations for these 
functions: 

'/'lOt = "00, 

"iot = 2uf , 

"OO* = "OOMlO, 

(f 2M = 2(u 00 - u 10 ip ), 
It is easy to find the solution of this system: 



"00 = ¥>20 - ¥>10- 



1 - 2u\t" 

f2Q = 0,2+ 2Ut, 



MOO 



(l-2ti?t)V2' 

</?io =0,1+1 u 00 dt, 
Jo 



V'o = [Wo«i - «o)(l - 2«5*) - "ov/ 1 - 2 < l 

£/ = Uq + «i(ai - 02). 

One can easily see that the function ipo(t) vanishes at the two points ti 
1/2mq and t* such that 



Obviously t* < ^ and the free singularities supports x = y> 10 and a; = </? 2 o 
merge at i = t* . In this case we have 



u 00 (t*) = Uq = U, u w (t*) 



U 2 



(u ^ 1 



u\ < 00. 



Thus in this example the mechanism of formation of a new shock wave consists 
not in turning over the inclined segment of the broken line, as in the preceding 
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example, but in the disappearance of this inclined segment due to increasing 
vertical segment. 

Subtracting the first equation from the third equation in (7), we obtain the 
following equation for the function ip: 

V>t - (no - 2V>ui)(l - 2Bi(A ¥ >/e)) 

or, denoting p = Aip/e = (tpo + V>o^i(r))/e, r = ipo/e, 

P 



ip' p = [uo - 2^ui](l - 2Si(p)), 



r 



1. 



Note that we can use the formula for ipo (and for the functions u no and u\o) 
only for t 6 [0, i* + S], where S > is any number such that S < t\ — t*. 

To obtain formulas that are global in t, we need to choose a number 5 and 
continue the functions uoo, uio, and ipo smoothly to the time t > t* + S so that 
the sign is preserved. Calculating the coefficient of m — 2i[>ui, one can easily 
see that p > for t < t*. Hence there exists a solution p — > p , where p is a 
root of the equation 

Bi(p) = 1/2- 

Let us consider the system of equations for the functions uq and U\. By the 
change uq = Uq-^/ui/uI, this system can be reduced to the single equation for 
ur- 

u u - 2u\{\ - 2Bi) = 0. 

Its solution has the form 



til (*,£•) = 



Ml 



l-2u\jl[l-2B l {p{r))]dt' 

Clearly, we have u\(t*,s) < u w (t*) (since J* (1 — 2Bi) dt' < t*). On the other 
hand, we have t > t* for p — > pq. Therefore, ipii T ) ~ * — 1 as r — > oo and hence 
(At/?) — > as r — > oo (i.e., for t > t* , e — > 0). This implies that for t > t* we 
have 

ui(t,e)=«io(t*) + o(l), e^O. 
We represent the above-constructed solution in the form 

u £ (x,t) = U6 St i(ipi - x) + (u - U)0 £t2 (<fi2 - x) 

+ [0e,l(Vl - X) - E ,2(f2 - X)](<fi2 - x)lH. 

Obviously, for t > t* , e — > 0, the first term approximates the shock wave 

u = U6{ai + Ut-x) 

and in this case the second term vanishes since tto — U — > and the third term 
vanishes since ip\ — <f2 = Ay — > 0. Recall that V'ot = u oo — 2?/;o u io f° r t < t* . 
In view of the equations, we can continue the function ip t for t > t* in the 
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form ipot = U . In this case the function uq — 2ipu\ is continuous uniformly in e 
for t = t* and we can show that the function p is a solution of boundary value 
problem 

P 



p=(l-2S 1 (p)), 



T 



1. 



The system of equations determining the weak asymptotic solution in this case 
also splits into separate equations. 

Now we briefly consider the problem of decay of nonstable shock waves. 

One can easily see that by setting vt, s {x, t) = u e (x, T—t), T > t* , we obtain a 
T-dependent family of weak asymptotic solutions of the equation vt — (v 2 ) x = 0. 

For t — the solutions of this family are shock waves (unstable for this new 
equation). The weak limit of these solutions for < t < T — t* is a shock wave, 
and for t > T — t* is a broken line consisting of two moving weak discontinuities 
into which the unstable shock wave splits at time t* (T) —T — t* (which is not 
unique). 

(c) Interaction of shock waves in the multidimensional case. Let us consider 
the two-dimensional nonlinear equation arising in the reservoir problem 

. — du , du 2 A du 2 
L[u] = — + A 1? — + A 2 — = 0. 
at ox i 0x2 

The above approach can be easily generalized to the case of an arbitrary 
dimension if the codimension of the front of the nonlinear wave is = 1. We 
assume that A\ ^ A 2 are positive constants. 

We choose the initial conditions as 

u \t=o = u o + u iS{t + ipi(x)) + u 2 0(t + ipi(x)), 

where x — (xi, x 2 ) 1 Ui are positive constants, and tpi(x) are the desired functions. 
We write T° = {x,i/>i(x) = 0}. 

Clearly, the curves T° are given initial positions of the fronts of two shock 
waves whose sum is just the initial condition. In addition, we assume that the 
curves T° are transversal to the vector field {A, V), A = (Ai,A 2 ) and T 2 is 
the cross-section of the (trivial) fibration over whose fibers are straight lines 
parallel to the vector A. In addition, we assume that the motion from the points 
of T 2 to the points of r° is in the direction of the vector A. In this case the fact 
that Ui,u 2 are positive constants is a sufficient condition of stability. 

If the functions ipi(x) are known, then the curves (level surfaces) Y\{x,t + 
■ipi(x) = 0} determine the fronts of shock waves at time t. 

Acting as before (see 2 (a)), we substitute the approximation 

u e (x, t) = u + uiO e ,i(t + -0i(x,e)) + u 2 9 £ ^ 2 (t + ip 2 (x,e)) 
into the equation and calculate the weak asymptotics of L[u £ ]. We obtain 

L[u e ] = S r t \l + (A, Vipi)(u\ + 2uq + 2u 2 Bi(Aip / e)) 



+8rt 



1 + (A, Vip 2 )(2u Q +u 2 + 2 Ul B 2 {Aip/e)) + Oxv (e), 



Alp = "02 - tpl ■ 
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Here the functions B\ and B 2 are the same as in item 2 (a). Formulas for the 
weak asymptotics in the multidimensional case are carried out in the same way 
as in the one-dimensional case 

Roughly speaking, the (two-dimensional) integral becomes an iterated inte- 
gral over the surface T\ and over the normal to this surface. The asymptotics 
of the integral along the normal is calculated in the same way as in the one- 
dimensional case, see [1]. 

It follows from our assumptions that the inequality Aip < holds for suffi- 
ciently small positive t. Hence we have Aip/s <~ —00 and for small t we obtain 
the following equations describing the system of noninteracting fronts: 

l + (l,Wi>(«i + 2u ) =0, l + (i*,VV>2}(u2 + 2wi + 2uo) = 0. (8) 

Clearly, these equations determine the (limit) functions tpko if the curves 
on which they vanish are given. 

Dividing these equations by |V?/>i| and taking into account the fact that, in 
view of our formulas, the waves travel in the direction of decreasing functions 
tpko(x), we can rewrite the last system as 

VjU = (A, ni)(ui + 2u ), V$ = (A, n 2 ){u 2 + 2m + 2u ), 

where is the normal (at a point) to Y\, Vnf is the normal velocity of this 
point. Clearly, the velocities of the points of the curve Y 2 are larger than the 
velocities of the points of the curve Y\ along the trajectories of the field (A, V), 
but the distance between the curves T* along the trajectory depends, in general, 
on the point. 

Therefore, since the shape of these curves is rather arbitrary, there may be 
no complete confluence of these curves at their contact. A new shock wave 
with summary amplitude u\ + u 2 is generated at the points of contact. This 
shock wave travels with its new velocity, and the solution may be of a rather 
complicated structure. To describe this wave uniformly in time, we shall seek 
the solution of the system 

l + (A,Vip 1 )(u 1 + 2u + 2u 2 B 1 (Aip/e)) = 0, (9) 
1 + (A, Vip 2 )(u 2 + 2u + 2uiB 2 (AV>A0) = 

in the form 

ipk(x,e) = ipko(x) + 0o(a#ki(0o/e)> (10) 

where 4>o = ip 2 o(x) — i>io(x)- Note that, in view of our assumptions on the 
geometry, instead of the coordinates (xi,x 2 ), we can introduce the coordinates 
(s, £), where s are the coordinates on T° and £ is a parameter on the trajectories 
of the vector field (I,V). 

Hence we, in fact, "calculate the distance" between the curves Yf (i.e., the 
differences Aip, <fi ) along the trajectories of the field (A, V). 
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Preserving, instead of j|, the notation (A, V), substituting (9) into (8), and 
taking into account (10), we obtain 

2(1, V^io>«2Si + (A, V0 o )|:(t^ii)[C/o + 2u 2 (Bi - 1/2)] - 0, (11) 
-2 Ul (l, W20) + 2(1, Vip 20 ) Ul B 2 

+ V^o)^-(t^2i)[^ + 2«i(B 2 - 1/2)] = 0. 

Here U — u\ + u 2 + 2u , t = 4>o/e. 

The further is similar to that in the one-dimensional case. Its first stage is to 
obtain an equation for the function p = Aip/e = (<j>o/e)(l + V ; 2i(t) — ipu (t)) = f 
r(l + <Pi(t)). Next, we calculate the limits of the functions Bk(p) as r — > 00 
(after the interaction) and find equations for the limit functions ip£ as well as 
equations for Vfci( r )- 

Subtracting the first equation in (10) form the second one and carrying out 
several calculations, we obtain the desired equation for p: 

j l-B 2 (p)\ 2u 2 U 2 , 2 Ul U x 



U\ + u 2 



1. 



.C7b - 2u 2 (B 2 - 1/2) l/i + 2«i(B2-l/2) 

One can show that the right-hand side of this equation (that differs, as one 
can see, from that in the similar equation in the one-dimensional case) also has 
a single root p and B 2 (p ) = 1/2 (and hence Bi(p ) = 1/2). 

Hence it follows from Eq. (8) that, for the same values of s for which a point 
of the curve "outruns" the curve , we have 

l + {A,V^)(u 1 +u 2 + 2u ) = 0, k = 1,2. 

This implies ipi = i\) 2 and for a given s, after the confluence of the curves, a 
wave with summary amplitude u\ + u 2 travels in the direction of A. Thus, for 
a fixed s, the dynamics of interaction in the direction of A is similar to that in 
the one-dimensional case. 

Here we do not write out the equations for ipki- They can be obtained in 
the same way as the similar equations in the one-dimensional case. 



3 Conclusion 

The problem of shock wave interaction in the one-dimensional case is presented 
in [2, 3, 5] not only for the Hopf equations but also for equations with sufficiently 
general nonlinearity. The formulas from Sec. 1 are derived there in more detail. 

Similarly, in the multidimensional case one can easily generalize our con- 
struction to the case of more general nonlincaritics, variable coefficients and 
amplitudes. 

For reasons of space, here we do not consider the problem of constructing 
definitions of weak solutions. This problem is discussed in [1, 4, 5]. In particular, 
in [5] a definition of a weak solution is constructed for KdV type equations 
admitting the zero dispersion limit for soliton type solutions. 
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